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The picture of the "multiverse" arising in diverse cosmological scenarios involves transitions be- 
tween metastable vacuum states. It was pointed out by Krauss and Dent that the transition rates 
decrease at very late times, leading to a dependence of the transition probability between vacua on 
the age of each vacuum region. I investigate the implications of this non-Markovian, age-dependent 
decay on the global structure of the spacetime in landscape scenarios. I show that the fractal di- 
mension of the eternally inflating domain is precisely equal to 3, instead of being slightly below 3 
in scenarios with purely Markovian, age-independent decay. I develop a complete description of a 
non-Markovian landscape in terms of a nonlocal master equation. Using this description I demon- 
strate by an explicit calculation that, under some technical assumptions about the landscape, the 
probabilistic predictions of our position in the landscape are essentially unchanged, regardless of 
the measure used to extract these predictions. I briefly discuss the physical plausibility of realizing 
non-Markovian vacuum decay in cosmology in view of the possible decoherence of the metastable 
quantum state. 



I. INTRODUCTION 

According to today's accepted cosmological data, the 
universe is now undergoing accelerated expansion with 
an approximately constant Hubble rate H now . However, 
models of string theory suggest that this accelerating 
state may be merely a metastable vacuum that is des- 
tined, after a long time, to decay via quantum tunneling 
into other states with different values of H . The recently 
developed paradigm of "string theory landscape" [l| in- 
volves a very large number of metastable vacua, corre- 
sponding to local minima of an effective potential in field 
space. The value of the potential at each minimum de- 
termines the effective Hubble rate in the corresponding 
vacuum. A similar scenario combining inflationary evolu- 
tion and tunneling was proposed earlier in Ref. [3] under 
the name of "recycling universe." In all these scenarios, 
the universe becomes a "multiverse," that is, an infinite 
ensemble of large, causally disconnected spatial regions. 
Some of these regions contain galaxies and stars, while 
other regions are undergoing inflation and generating new 
vast domains of space. Each spatial domain may be in 
a metastable vacuum state with a sufficiently long decay 
time, so that reheating can occur and the standard cos- 
mological evolution can proceed before the transition to 
a different vacuum state. 

The theory allows in principle to determine the set of 
possible vacua but does not predict our position in the 
landscape with certainty. After many transitions, the 
position of our observable patch of the universe in the 
landscape becomes random. Nevertheless, one would like 
to explain the present value of the cosmological constant 
and possibly other observables. Therefore one attempts 
to calculate the probability of being in a vacuum of a given 
kind, for a "typical" observer. It is notoriously difficult 
to formulate an unambiguous and well-behaved measure 
on the set of all possible observers such that the "typical" 
observers are selected without bias; see e.g. d, H[ for a 



recent discussion and Refs. [U, [f| 0, [1] for reviews of the 
proposals of observer-based measure. 

In this paper I study a different aspect of the measure 
problem. All currently proposed measures are based on 
the assumption that the decay of a metastable state pro- 
ceeds independently of the individual age of that state. 
In other words, it is assumed that the random process 
of transitions between different states in the landscape is 
a Markov chain. Markovian transition probabilities are 
determined only by the current state and have no mem- 
ory of previous transitions. (The "memory" effect due to 
bubble collisions @ does not modify transition probabili- 
ties.) Vacuum decay proceeds through bubble nucleation 
and is normally described via the nucleation rate per unit 
4-volume US 111, 



=<3(l)ff 4 exp 



(1) 



where Si is the relevant instanton action and H is the 
Hubble rate of the parent vacuum (I use the Planck units 
throughout the paper) . The transition rates between dif- 
ferent metastable vacua can be considered (in principle) 
known in a given model of the landscape. For a fixed 3- 
volume V, the probability of nucleating no bubbles after 
time t is exponentially small, <x exp[— T^ D ^ Vt\. 

A statistical description of evolution in the landscape 
can be obtained @, [l2j] by considering the fraction f a (t) 
of the comoving volume occupied by bubbles of type a 
at time t. One can approximate the transition to a dif- 
ferent vacuum as a series of random nucleation events, 
each event resulting in an instantaneous conversion of a 
volume H~ 3 of vacuum type a into the same volume of 
vacuum type (3. The rate of this conversion per unit time, 
denoted n a —,/3, can be computed according to Eq. ([T]) 
with appropriate normalization factor, 



0(l)H a exp 



(2) 
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Defining for convenience n a ^ a = 0, one then writes the 
master equation describing the evolution of f a (t) , 

-jT = ( K P~>afp ~ K-a^pfa) ■ (3) 

at p 

This equation can be solved with initial conditions f a (0) . 
Some measure prescriptions based on the comoving dis- 
tribution were proposed in Refs. [HI, [I3> EH • 

Another useful distribution is the 3- volume V a (t) of 
spatial regions within bubbles of type a at time t. The 
evolution equation for V a (t) differs from Eq. ^ by the 
volume expansion factors, 

= 3H a V a + ^ («/3-«*//3 - K a ^pf a ) (4) 

at p 

= J2 M ^Vf3, (5) 



where the matrix M a p is defined by 

M a/3 = ^3H a - ^ K a ^\j S a p + Kp^a. (6) 

The volume-weighted master equations are used in 
volume-based measure prescriptions (e.g. Refs. [l2l fl6|). 

All the existing measure prescriptions depend on the 
properties of the late-time behavior of the distributions 
f a (t) and V a (t). The late-time asymptotics of the so- 
lutions of Eqs. ©-([I]) are always exponential. For in- 
stance, the volume distribution has the late-time asymp- 
totics V a oc c a e l1 , where 7 > is the dominant eigen- 
value of the matrix M a p. The values of the coefficients 
c a are determined by the right eigenvector of M a p cor- 
responding to the eigenvalue 7. 

Recently, Krauss and Dent [TtJ called attention to 
the fact that the decay of metastable states becomes 
subexponential at very late times. In typical quantum- 
mechanical metastable systems in d-dimensional space, 
the probability of not decaying (the "survival probabil- 
ity") initially decreases exponentially as e _rt , where T 
is the decay rate, but eventually starts falling off as t~ d 
after a (very long) crossover time T ~ 5r _1 In (E/T), 
where E is the energy difference between the metastable 
state and the final stable state. Effectively, the tunnel- 
ing rate for all transitions between states goes to zero as 
T(t) oc i _1 after a (state-dependent) crossover time. It is 
important to note that the transition dynamics depends 
on the "age" of the current state, i.e. on the time elapsed 
since the last transition. With this modification, the 
transition process becomes a non-Markov random walk, 
and Eqs. ©-(Ill) no longer apply. In particular, the late- 
time asymptotics of the bubble distributions fpit) and 
Vp(t) are no longer purely exponential. For this reason it 
is interesting to investigate the implications of the non- 
Markov transitions for the measure calculations, which 
depend in an essential way on the late-time behavior of 
fpit) and Vp(t). 



In this paper I study the evolution of the landscape 
assuming that the late-time asymptotic of the survival 
probability becomes subexponential at a state-dependent 
crossover time. The main results of this first study are as 
follows. I show that the fractal dimension of the inflating 
domain is exactly equal to 3, while it is always slightly 
below 3 in Markovian models. Then I develop an explicit 
non-Markovian description of the transition dynamics in 
terms of a master equation that is nonlocal in time. Using 
that equation, I derive the late-time asymptotics of the 
volume distributions V a (t) using the proper time coordi- 
nate t. The results show explicitly, within a controlled 
approximation, that the volume ratios V a (t)/Vp(t) ap- 
proach a constant at late times and are approximately the 
same as those computed within the Markovian situation, 
except for the volume in bubbles of type having the 
largest Hubble rate Hq = max a H a . The bubbles of type 
now entirely dominate the volume of the universe at a 
fixed time t, whereas their volume fraction was large but 
finite in Markovian scenarios. These results (obtained us- 
ing the proper time gauge) applied to landscapes where 
a single vacuum type has the largest Hubble rate of all 
available vacuum types. I also show that the comoving 
volume distributions remain essentially unchanged in the 
non-Markovian regime. This suggests that the results 
obtained in any measure prescription (whether volume- 
based or worldline-based) do not need any modification 
in view of the modified late-time decay. I conclude with 
a brief discussion of the viability of the non-Markovian 
assumption in the cosmological context. 



II. NON-MARKOVIAN SIERPINSKI CARPET 

I begin by examining the global structure of the space- 
time undergoing non-Markovian vacuum decay. A par- 
ticular version of the random Sierpihski carp et, or "in- 
flation in a box," was considered in Ref. [18| as a dras- 
tically simplified toy model mimicking the global geom- 
etry of such a spacetime. In this model, time elapses 
in discrete steps, and the space is reduced to a two- 
dimensional square domain < x, y < 1, where x,y are 
the comoving coordinates. The entire initial Hubble-size 
domain is assumed to be initially inflating. To imitate 
inflation during one time step, one subdivides the ini- 
tial inflating square into N x A equal sub-squares of size 
A -1 x A -1 ; at the next step, each sub-square will again 
have the Hubble proper size. Then one randomly marks 
some of the smaller squares as "thermalized," assuming 
that each Hubble-size inflating square continues inflation 
with a probability q (where < q < 1) and thermal- 
izes with probability 1 — q, independently of all other 
squares. The selection of thermalized squares concludes 
the simulation for one timestep. At the next timestep, 
the same procedure of subdivision and random thermal- 
ization is applied to each Hubble-sized inflating square, 
while the "thermalized" squares do not evolve any fur- 
ther (see Fig. [!}. This process is continued indefinitely 
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Figure 1: First steps in the construction of a random Sier- 
piriski carpet with N = 5 and q = 5/6. 

and generates a fractal set of measure zero consisting of 
points that never enter an y th ermalized squares (called 
the "eternal points" in Ref. [19| where rigorous definitions 
are given). This set represents the eternally inflating sub- 
domain of the spacetime. Under the condition N 2 q > 1, 
the fractal dimension of the eternally inflating domain is 
7 = 2 + In g/ In TV > 0, and future-eternal inflation occurs 
with a nonzero probability [l8j . 

As formulated, the model is Markovian since the ther- 
malization probability at each step is independent of the 
age of the inflating square. The probability of remain- 
ing in the inflationary regime (the "survival probability") 
after t time steps is q t — e~ at , where a = ln|. Let us 
now modify this toy model by assuming that the survival 
probability is given by a function S(t) that interpolates 
between the initially exponential falloff S(t) — e~ at for 
t <C T and the power-law asymptotic S(t) « Sot~ p for 
f > T, where p is a fixed constant and T is the crossover 
time. We would like to compute the fractal dimension of 
the set of eternal points in this non-Markovian model. 

Let us denote by X(t) the probability of the presence 
of at least one eternal point within an inflating square 
at time t. The quantity X(t) can be computed explic- 
itly, but it is sufficient for the present purposes to obtain 
the asymptotic value of X(t) at t —> oo. Since the ther- 
malization probability per step goes to zero at late times, 
the value of X(t) approaches 1 as t — > oo. More precisely, 
X(t) is the nonzero solution of the equation 

l-X(t)=p(t) + (l-p(t))(l-X) N \ (7) 

An approximate solution of this equation for pit) <C 1 
is f — X w pit). Since p(t) — > as t — > oo, we have 
X(t) — > 1. Hence the average number of inflating squares 
containing at least one eternal point at a late time t is 
~ S(t)N 2t , while the linear size of each square is A -1 . 
So the fractal dimension of the eternally inflating set is 

7 = 2- lim H = 2+ lim — y — '- = 2. (8) 

It can be shown that the eternally inflating domain con- 
sists of an infinite merged cluster when it is formed as a 
random Sicrpinski carpet with fractal dimension 2. It is 
important to note that the eternal set still has measure 
zero because every comoving point will reach thermaliza- 
tion with probability 1. 



By analogy, one can investigate the eternally inflat- 
ing domain in a three-dimensional space and conclude 
that its fractal dimension is 3. A quick argument leading 
to this conclusion consists of estimating the growth of 
the 3-volume of the inflating domain as V(t) oc e 3Ht t~ p 
in the regime of power-law decay at late times. A do- 
main growing as V(t) oc e lHt is interpreted as a lacunary 
fractal with dimension 7 [l8|, [2(|, regardless of subex- 
ponential corrections. Therefore, the fractal dimension 
of the inflating domain is always equal to 3 in the non- 
Markovian case. This is only a small correction to the 
results obtained in typical scenarios of eternal inflation 
where the fractal dimension is very slightly below 3 (see, 
for instance, Refs. [13, EH). Therefore, the global geome- 
try of the spacetime is not significantly modified in these 
scenarios even if the late-time decay is subexponential. 



III. EVOLUTION IN A NON-MARKOVIAN 
LANDSCAPE 

The next issue is whether the results of applying the 
various landscape measure proposals are modified when 
non-Markovian decay is assumed. In this section I derive 
the suitably modified versions of Eqs. © - ® and ob- 
tain their late-time asymptotics. Since all the different 
measure proposals require computing the late-time be- 
havior of these same evolution equations, the results of 
the present calculation will be equally relevant to every 
measure proposal. 

To describe the evolution of spacetime in a landscape 
scenario with non-Markovian transitions, one needs to 
specify the transition rate T a —,p (t) between vacua a and 
P as a function of the age t of the parent vacuum a. 
The precise form of T a —,p(t) will be model-dependent 
except for the properties T a —>p(t) w n a ^fj = const 
for t < T a ^f3, where T a ^p is the crossover time, and 
^a~>ii{t) — * for t ^> T a ^p. For simplicity I will as- 
sume below that the crossover time T a —,p = T is inde- 
pendent of a and (3. Without this technical assumption, 
the analysis will be more complicated without yielding 
significantly different results. If transitions a — > (3 have 
different crossover times T a ^p, the results of the present 
analysis will be approximately applicable at sufficiently 
late times t such that t 3> T = max Qii a T a ^p. 

Since transition probabilities depend on the age, it is 
not sufficient to consider the probability distributions 
fp[t) and Vfs(t) mentioned above. One needs to intro- 
duce more detailed distributions that include information 
about the times of the previous transitions. 



A. Volume distributions 

I first consider the volume distribution. Assume for 
convenience that there is a single initial bubble of type 
ao formed at time t = with unit volume, and that we 
are interested in describing only the evolution of the inte- 
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rior of the initial bubble and any bubbles nucleated in it. 
(The case of several initial bubbles is a straightforward 
extension.) Let V a (to,t)dto denote the volume at time t 
of bubbles of type a that were formed at an earlier time 
between to and to+dto- By definition, we set V a (to, t) = 
for to > t. The volume remaining from the initial bubble 
could be included in V ao (to,t) as a contribution of the 

form 6(to)Vdl (t), but it is techically more convenient to 
exclude the initial bubble from V a (to,t) and to account 

for its volume Va„\t) separately. The quantity Va® (t) 
represents the proper volume that remains from the ini- 
tial bubble and has not decayed by time i. The volume 
Vao\t) of the initial bubble grows with the rate 3H ao 
and decreases due to nucleation of other bubbles: 

= 3H ao V£Ht) -J2r ao ^(t)vM(t). (9) 

Integrating Eq. §§§ with the initial condition vi° ' (0) = 1 , 
we find 



V£\t) =exp[3H ao t]S ao (t), 



S a (t) = exp 



5> 

p 



>0{t')dt' 



(10) 



(11) 



The auxiliary function S a (t) is the survival probability 
of a bubble of type a and age t. 

The evolution equation for V a (to,i) accounts for the 
growth of volume at rate 3iJ Q , age-dependent decay into 
bubbles of different kinds, and age-dependent nucleation 
of zero-age bubbles of kind a from other bubbles (includ- 
ing the original bubble): 

dV a (t ,t) 



3H a V a (t , t) 



J2 T cc^p(t-t )V a (t ,t) 

P 



5{t-t ) / dt J2 



Of 



P 

+ 6(t-t )T ao ^ a (t)vW(t). (12) 

The factors 5(i — to) account for the fact that bubbles nu- 
cleated at time t have zero age at that time and therefore 
contribute to the distribution V a (to, t) only at to = t. 

Noting that V a (to,t) with t ^ to is decoupled from 
other Vp(to,t), we have (for t > to) 

V a (t ,t)=V a (to,to)exp[3H a (t-t )}S a (t-t ). (13) 

It remains to determine the function V a (to, to) = U a (to). 
Integrating Eq. (|12[) in t over an infinitesimal interval 
around t = to and using Eq. (|13p and the condition 
V a (t 07 t) = for t < to, we obtain a closed system of 
integral equations for U a (t), 

M<) =Y] fdtoU^e^-^Spit - t )T p ^ a (t - t ) 

(14) 

+ T ao _ tQ (t)e 3H ~o t S ao (t). (15) 



It remains to determine the asymptotic behavior of the 
functions U a (t). 



1. Markovian regime 

I first consider times t before the crossover time scale, 
< t < T. At these times, the behavior of the system is 
(approximately) Markovian, and one expects to recover 
the standard equations (j4|). In Eq. (| 15(1 we may approx- 
imately set 

T/3^, a {t -t ) w Kp->a = const, (16) 
Spit -to) « exp [- (i - to) I>] , (17) 

where we denoted by Tp = J2 a K P^a the total decay rate 
of the vacuum type (3 in the Markovian regime. Then 
Eq. (fT5]) is rewritten as 



U a (t) 



. r 

J- rv n — * <y 



dtoUp{to)e 



,(3H ao -r ao )t 



(18) 



Although this system of equations appears to be nonlo- 
cal in time, it can be reduced explicitly to a Markovian 
system. We pass to new variables 

V a (t) = fdioU a (io)e^- T ^-^ +S aao e^o-^o)K 
Jo 

The quantity V a (t) represents the total volume inside 
bubbles of type a at time t integrated over the bubble 
ages and also including the volume of the initial bubble. 
The variables U a {t) are expressed through V a (t) as 



U a (t)=e^- r ^d t 



~(3H a 



= V a -(3H a -T a )V a . (20) 
Hence the volumes V a (t) satisfy the differential equation 

V a = (3H a - T a ) V a +J2 Kp^ocVp = M ^ V f J ( 21 ) 
P P 



with the initial condition ^(O) 



which is equiv- 



alent to Eqs. (U]); the matrix M a p is defined by Eq. 
The late-time asymptotic of solutions is exponential, 



V a (t) 



alt 



(22) 
It 



where 7 is the largest eigenvalue of the matrix M a p 
is important to note that 7 > 3Hp — Tp for all /3. 

The eigenvalue 7 and the corresponding eigenvectors 
of M a p can be estimated explicitly under some technical 
assumptions. To be specific, let us denote by Ho and Hi 
the first and the second largest values among all the H a , 
and let us assume that the nucleation rates are small, 



K a ->p <C Ho — Hi for all a, j3. 



(23) 
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Since the nucleation rates are typically exponentially 
small, one can disregard terms of higher order in K a ^. 
Then the matrix M a p can be represented as a diagonal 
matrix S a p (3H a — T a ) with a small perturbation of order 
K a ^p , and the dominant eigenvalue is found by standard 
perturbation theory as 



7 



3H -T + Y: rr° K °r + O(^ Y (24) 



Q#0 



3Hq — SH a 



The second-order term in 7 will play a role below. 

The coefficients c a in Eq. ([22| are proportional to the 
components of the (right) dominant eigenvector r a0 of 
M a p, so that c a /c/3 — r a Q/rpQ. The ratios of components 
of the eigenvector r a o can be found approximately as 



KQ- 



r a o 



r 00 3H - 3H C 



+ 0(K 2 a ^), a^O. (25) 



It is useful to compute also the absolute normalization 
of the coefficients c a , which will yield an explicit late- 
time asymptotic V a (t) = c a e 7 * as a function of the initial 
conditions V Q (0) = 8 aao . The time-dependent solution 
V a (t) can be decomposed as 



V a {t) = v n r an e lnt , 



(26) 



Hence, a suitable normalization of the eigenvectors is 
roo = 1, r a0 = nTT K °^* TT + 0(K%_>p), a^O; 



'00 — It ' 



aO 



3Hq — 3H, 



(31) 
I. 

(32) 



Now we may compute explicitly 

c a — v r a0 — r a0^a 0- (33) 

The full solution U a (t) can be written as 

U a (t) = ^l aon r a n (in - 3H a + T a ) e 7 "*, t < T. 

n 

(34) 

The late-time (but still Markovian) behavior of U a (t) is 

U a (t) « r a0 l aoQ (7 - 3H a + T a ) e 7 *. (35) 

Although the absolute values of the coefficients c a de- 
pend on the initial conditions, the ratios c a /cfj do not. 
This is the standard property of Markovian models: the 
late-time asymptotics do not depend on the initial con- 
ditions. 



where 70,71, ... and r aQ , r Ql , ... are the eigenvalues and 
the corresponding (right) eigenvectors of M a p. The late- 
time behavior of V a (t) is dominated by e 70 ' , where 70 = 7 
is the largest eigenvalue. 

The coefficients v n are found by decomposing the ini- 
tial condition vector V a (0) in the basis {r Qn }, 



(27) 



The coefficients v n are computed as the products of the 
left eigenvectors l an , n — 0, 1, ... of M a p with the initial 
condition vector V a (0) = 5 aao , 



5^UV r a (0)=/ o 



(28) 



where we assumed that the dual bases {l an } and {r an } 
are normalized, 



El r = S 
L am 1 an u n 



(29) 



We are interested only in the coefficients V a o correspond- 
ing to the dominant eigenvalue 7 = 70, so vq — l ao o- 
The vector l a Q is determined perturbatively under the 
assumption (T2"3"l) through the ratios 



laO _ K a^0 

loo 3i?o ~ 3H a 



+ 0(4^), «^0. (30) 



2. N on- Markovian regime 

Having computed the early-time behavior of U a (t), I 
now consider the asymptotics of U a (t) at late times t for 
which the survival probabilities S a (t) are subcxponcntial. 
Since the decay rate is the logarithmic derivative of the 
survival probability, it follows that T a ^p(t) oc t^ 1 at 
those times. To simplify calculations, I assume that 



S a (t)T a ^i3 = R a (t)K a ^.p for all /3, 



(36) 



where the function R a (t) describes the transition from 
the Markovian to the non-Markovian regime as 



Ka{t) -\e^r a T}(^, t>T, 



(37) 



where T is the crossover time and p a > are constants of 
order 1. (For the cited examples of subexponential decay 
with S a (t) oc t~ 3 one will have to set p a = 4.) The as- 
sumption of a common time profile R a (t) and a common 
crossover time T, independent of the vacuum type a and 
of the decay channel a — > (3, may be insufficiently precise 
in some scenarios. Here I employ this technical assump- 
tion as a first step towards a more complete calculation. 

Let us first determine the ansatz for the asymptotics of 
U a {t) by examining Eq. (fT5|) . Since U a (t) receives con- 
tributions from all the subexponentially decaying states 
(3 a according to Eq. (fTS"]) . the late-time asymptotics 
of U a (t) must grow at least as fast as the fastest-growing 
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function among e 3Hat R a (t) for all a. Hence, the expo- 
nential part of the asymptotic is U a (t) oc e 7 *, where 7 is 
not less than 3Ho and Hq is the largest available value 
among H a . However, the function U a (t) cannot grow 
faster than e 3Hot , i.e. as e 7 * with 7 > 3Hq, because in 
that case the integral in line (fl4|) is dominated by to ~ t 
(recently nucleated bubbles) where the survival probabil- 
ities Sp(t — to) are Markovian. So the Markovian calcula- 
tion leading to Eq. (fl8|) still holds and yields the contra- 
dictory result 7 = 7 ~ 3H ~Tq < 3H - Hence, 7 = 3H - 
We need to allow for the possibility that U a (t) contains 
also a subexponential asymptotic, U a (t) oc e 3Ht)t Q(t), 
where Q(t) is a subexponential function decaying not 
faster than R a (t) at t > T. (Below I will show that 
Q(t) oc Ro(t), but at this point the behavior of Q(t) is 
not yet determined.) Thus, the late-time asymptotics of 
U a (t) are of the form 



U a (t) ~q a e 3Hat Q{t), t >T, 



(38) 



while the Markovian behavior was determined above in 
Eq. ([34]) . The task at hand is to determine the coeffi- 
cients q a and the function Q(t) for the non-Markovian 
asymptotics (|3"5j) . 

Let us define the auxiliary quantities 



integrals are estimated as follows. The first integral, 

i-T 



dtoUoi^e^-^Roit-to), 



(42) 



is dominated by the contribution of to ~ because L/o(fo) 
in the Markovian regime grows as e 7 * , while 7 < 3Hq. 
Using Eq. ([51]). we find 



dtoUo(t )e 3Ha{t - to) R (t - t a ) 
e^R(t)^v n r J\^ Ho + To . (43) 

~~ JflO — In 



The sum in the last line can be estimated without ac- 
tually computing all the eigenvectors ro« by noting that 
3Hq — 7„ 3> T for n ^ 0, and thus the factor 



in - 3H + r n 



3H Q - 7„ 
For n — this factor is negligible, 



-1 + 0(k q ^), 71 ^ 0. (44) 



7 - 3H + r ^h-x v fAK , 
= O{H K a ^p), (45) 



W a (t) = / dtoUMe^'-^R^t-to) 
Jo 

+ S aao e 3H « t R a (t), (39) 



so that Eq. (fTS")) becomes 



m*)=E 



We will first determine the asymptotics of the quantities 
W Q (t) for t;>T. 

The definition of W a (t) involves an integral over to 
that needs to be estimated. It is convenient to estimate 
it separately for a^O and a = 0. For a/0, the function 
U a {to) grows as e 7 * until t — T; subsequently U a (to) 
grows even faster, as e 3Hot . This function is multiplied by 
a decay factor e~ 3Hat ° R a (t ~ to) that never compensates 
the growth of U a (to) if a ^ because 7 — 3H a ;» T a 
for a / 0. Therefore, the integral over to is dominated 
by the contribution near the upper limit to ~ t where 
R a (t— to) is Markovian while U a (t) oc e 3Hat . One obtains 
the asymptotic estimate 



W a (t) 



q a Q(t) 



3Hq — 3H a + r c 



a + 



(41) 



where the term oc e 3Hat R a (t) can be disregarded since it 
is exponentially smaller at late times. 

Estimating the quantity Wo(t) requires somewhat 
more work. One needs to split the integral in the defini- 
tion of W a (t) into three subintervals [0, T] , [T, t — T], and 
[t — T, t] where different factors in the integrand have ei- 
ther Markovian or non-Markovian behavior. These three 



where we used Eq. (jM)) . By splitting off the n = term 
from the sum in Eq. (|43|1 . one now obtains 



E7„ - 3H a + T , 
v n r 0n — 577 «-> v n r „. (46) 
3H - 7„ ^ 



(40) The last sum can be evaluated using Eq. (|27|) . 



Vo(0) = v n r 0n = v Q r o + v n ro n = So ao , (47) 

n n^O 



and we find 



E 



7„ - 3H Q + To 



v -Soa - (48) 



3H - 7„ 

Hence, the expression (|4"3"|) is estimated as 

e 3ffot i?o(t)(i;o-^o Q o)- (49) 
The integral over the second interval, 



t-T 



dtoUo{to)e 3Ho(t - to) Ro{t-to), (50) 



involves both Uoito) and Ro(t — to) in the non-Markovian 
regime. We find 



t-T 



d£ EMio)e 3ffo( *~ to) .Ro(t-io) 



9o e 



3iY t 



t-T 



dt Q(t )Ro(t-t ). (51) 
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Since both Ro(t) and Q(t) are decaying functions, we 
may estimate the integral in Eq. (|5ip as the sum of the 
contributions from intervals of order T at the two ends 
to = T and to = t — T, 

q e 3Hot [R (t)Q(T) 0(T) + Q(t)R (T) 0(T)} . (52) 

This precision is sufficient since these terms will not play 
a significant role in the final result. 

The integral over the third interval involves the Marko- 
vian Ro(t — to) and is dominated by to ~ t, 



dtoU (to)e 3H ^-^Ro(t~to) » ^-Q(t)e 3H °\ (53) 

i-T 1 



where we disregarded e" r ° T <C 1. (Note that T T > 1.) 

Putting together the contributions of the three inter- 
vals as well as the last term in Eq. (f3"9")l . we obtain 



W (t) « e 3Hot R (t)v + qoe 3Hot T^Q{t) 



+ qoe 3H ° t R (t)Q(T)0(T), 



(54) 



where we disregarded 

Q(t)R (T) 0{T) < qoe 3Hat T^Q{t) (55) 

because 

Ro{T)0{ToT) = e- r " T 0(T Q T) < 1. (56) 



Finally, we substitute the ansatz (f3"B")) and the esti- 
mates ([54]) into Eqs. (gDJl for f7 a (i). In the limit 
< 3> T, we may divide through by the factor e 3Hot Q(t) 
and obtain a system of equations for q a and Q(t), 



q a = k/3^ q 

/3 



3i/ -3^ + r /3 ' 

3F - 3H/3 + T 



(57) 



+ K ^« N + g Q(T)O(T)] lim 



Ro(t) 

Q(t) • 



(58) 



This is an inhomogeneous linear system for {q a }- 

Let us consider the possible values of 
limt^oo Ro(t)/Q(t) that show whether Q(t) is asymp- 
totically dominant over Ro(t). Since Q(t) in any case 
does not decay faster than Ro(t), there are only two 
possibilities: either the limit is zero or it is nonzero. I 
will now show that this limit must be nonzero. 
If ]im t -Kx>Ro(t)/Q(t) = 0, we rewrite Eqs. (|S"7)- 

as 



E K f3-^a; 
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3H Q -3H P + T^ 
Passing to auxiliary variables 

q a 



a = 0, 1, 



3Hq — 3H a + T a 



(59) 



(60) 



we find 



3Hqs c 







(61) 



Since the largest eigenvalue of M a p is 7 < 3Ho, it follows 
that 3Hq is not an eigenvalue of M a/ g. Hence, the only 
solution of the homogeneous system (f59f is q a =0. This 
contradicts the assumption that q a e 3Hot Q(t) is the lead- 
ing asymptotic of U a (t). Therefore, Q(t) decays exactly 
as Ro(t) at late times. 

Since Eq. l[5"8|) depends only on the ratio Q(T)/Q(t), 
the normalization of the Q(t) could then be adjusted such 
that limt^oo R {t)/Q(t) = 1. The value Q(T) is of order 
e -r T ^ ue ^ ^ e continuity requirement 

3ffoT nm _ „ e 3 7 T 



U a {T) ~ q a e Q(T) ~ c a e 



(62) 



Therefore, the term q Q(T)O(T) in Eq. ((SSJ) is exponen- 
tially small and can be neglected. We note, however, that 
its magnitude depends on the initial conditions through 
the coefficient c a ~ O(T ), which introduces, strictly 
speaking, an exponentially small dependence on initial 
conditions, of order 0(ToT)e~ raT . 

Finally, we rewrite Eqs. (j5Tj) - (fB"8")) through the vari- 
ables s a as 



3H S a - ^2 M af3S(3 = V K ^c 




(63) 



This is an inhomogeneous system of equations with a 
nondegenerate matrix, and so the solution is unique. It 
follows that all s a are of order voK a ^/3, so an approximate 
expression for the solution is readily found as 



so 



v 



E 



K 0^O K O^0 



u 

VoKo^u 



'3Hq — 3H a 



3H - 3H 



(64) 
(65) 



The corresponding values of q a (neglecting higher orders 
of k„^) are 



qo « 



E 



3H — 3Hp 



q a w v Ko->a, a^0. 



(66) 
(67) 



We note that the solution depends on the initial bubble 
through i>o only in the overall normalization; the ratios 
qa/qp are independent of vq- 

Having determined the auxiliary quantities U a (t), we 
can now compute the non-Markovian volume distribution 
V a (t) as 



V a (t)= / dtoU a {to)e 3H ^^S a {t-to) 
Jo 



„3H an t 



ot s ao (t). 



(68) 



For a^0, the integral in Eq. (f68|) is dominated by to ~ 
t, which yields a term cx e 3Hat , so the second term in 
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Eq. (|6"5)) is negligible. Hence, by setting Q(t) « Ro(t) 
and q a = UrjKo-»a one obtains the estimate 



V a (t) 



V K ^ a 



3-ffo — &H a + Fq, 



Ro(t)e 



3H ( ,t 



a^O. (69) 



We note that the ratios of volumes V a (t)/Vp(t) are in- 
dependent of the initial condition parameter w and of 
time, indicating a "stationarity" of the solutions V a (t) 
with a / 0. Moreover, these ratios are equal to the ra- 
tios obtained in the Markovian regime, 



km ^4 



a,j9 ^ 0. 



(70) 



Tke imprecision in the above equality is exponentially 
small, of order 0(ToT)e~ r ° T , as noted before. (To sim- 
plify calculations, we also carried an imprecision of order 
Ka-tp/iHo — Hi) in the expressions for c a , but Eq. (|TD|) 
also carries that imprecision. This limitation is due to 
the approximations adopted in the present paper.) 

It remains to compute Vo(t). For a = 0, the estimation 
of the integral in Eq. (|68[) proceeds similarly to the argu- 
ment leading to Eq. (|54|) . except that Ro(t) is replaced 
by So(t) which decays slower. The result is 



V (t) ' 
Since at large t 



o 3H t 



[v S (t) + q r^Q(t)] . (71) 



Q(t) « R (t) = So(t) r °^ a{t) « So(t), 



K - 



the dominant asymptotic for Vo(t) for t 3> T is 
V (t) « e 3ffot f 5 (t). 



(72) 



(73) 



other types of bubbles less likely in very old regions. The 
absolute dominance of bubbles of type will set in after 
time T. This is different from the Markovian situation 1 
where bubbles of type dominate with a finite (but very 
large) ratio, 



l im \ * 

t ^oo ^/Markov^) 



£0 

Cq 



— 3iJ Q 



> 1. 



(75) 



Thus the qualitative picture of the distribution of volume 
in space has changed due to the non-Markovian decay, 
but the change is not drastic. This conclusion is similar 
in spirit to that obtained in Sec. [TTl where the fractal 
dimension of the eternally inflating domain was modified 
from 3 — e, where e <C 1, to exactly 3. 

On the other hand, the 3- volumes V a (t) within other 
types of bubbles a ^ grow proportionally to each other, 
and the ratios V a /Vp are almost the same (up to exponen- 
tially small corrections) as those obtained in a Markovian 
calculation. Therefore, any measure prescription that de- 
pends on the asymptotic ratios of volumes, V a /Vp, will 
give unchanged predictions as long as one asks about the 
volumes of bubbles of subdominant types (a ^ 0). Since 
the bubbles of type (presumably, with a Planck-scale 
Hq) are not especially interesting observationally, one can 
conclude that a possible non-Markovian decay has no ef- 
fect on predictions obtained via any measure prescrip- 
tions based on volume ratios. 

The considerations in the present paper are limited to 
proper time gauge and to landscape scenarios satisfying 
the assumptions (|23|) . The methods developed here are 
applicable to landscapes of any type, and future work will 
show whether the conclusions hold in more general cases. 



C. Comoving distributions 



B. Discussion 



We will now interpret the results of the calculation in 
the previous section. Since Q(t) <C So(t) at late times, 
the volume Vq (t) within bubbles of type grows asymp- 
totically faster than all other V a (t) for a ^ 0, 



lim 



Vo(t) 

— oc lim 

V a (t) t^oo 



S (t) 
R (t) 



= oo, a^0. (74) 



This indicates that the 3- volume at time t is entirely dom- 
inated by the bubbles of type 0, which we have labeled 
as those having the largest Hubble rate Hq = max a H a . 
Moreover, since the integral in Eq. (|6"5|) for a — is dom- 
inated by to s=s 0, it follows that almost all of the volume 
in bubbles of type at time t is in the very old regions 
of type 0. These regions of type either belong to the 
original bubble (if a = 0), or were nucleated early on 
(if ao ^ 0) and, by chance, have remained without de- 
cay for almost all of the time t. This dominance does 
not depend on the initial conditions and is due to the 
fact that non-exponential decay makes the nucleation of 



I now turn to considering the comoving distribution. 
One can define the distribution f a (to,t)dto as the frac- 
tion of comoving volume at time t in bubbles of type 
a that were formed at an earlier time between £o and 
t + dt . As before, we set f a (t 0l t) — for to > t; the 
volume remaining from the initial bubble is not included 
m fa{to,t) but accounted for separately as the function 
faj(t). The formalism and the equations for the dis- 
tribution / a (to,t) are quite similar to those developed 
above for the volume distribution V a (t , t) except for the 
absence of the volume growth factors H a . 



1 The 3-volume is not a gauge-invariant quantity, and statements 
about dominance of 3-volume at fixed time depend sensitively 
on the choice of the time variable H8|| . In particular, in Marko- 
vian models the 3-volume is not dominated by fastest-expanding 
bubbles if one chooses the e-folding time r = In a as the time 
coordinate. A similar gauge dependence is expected in the non- 
Markovian case. The present calculation focuses on the effects 
of non-Markovian decay, which are arguably more pronounced in 
the proper time gauge. 



9 



Instead of writing out the equations and the solutions 
f° r fa (to 7 1) , a simple consideration suffices to show that 
non-Markovian effects are irrelevant for the distributions 
of comoving volume. The comoving volume fractions 
f a {t), defined regardless of age, exponentially quickly 
become constant because the total comoving volume is 
conserved, and the dominant eigenvalue of the relevant 
Markovian matrix is equal to zero. The nucleation of 
bubbles will be always dominated by new bubbles rather 
than by "aged" comoving volume, simply because the co- 
moving fraction of the aged volume quickly goes to zero. 
In fact, the "aged" comoving volume has a smaller nucle- 
ation rate, r a _ >| g(t) — > for t — > oo, and therefore plays 
an even less significant role in nucleation of new bub- 
bles as in Markovian models. This is in contrast to the 
situation with the volume-weighted distributions, where 
the aged volume is rewarded by an exponentially large 
extra growth factor e 3Hot compared with the new vol- 
ume that grows slower, as e^ 3H °~ r °' t . Therefore, the 
non-Markovian decay law will introduce only a vanish- 
ingly small correction to the predictions obtained through 
comoving- volume measure prescriptions. 



IV. IS AGE-DEPENDENT DECAY 
COSMOLOGICALLY RELEVANT? 



A subexponential asymptotic at late times is a generic 
feature of quantum-mechanical systems. This feature 
can be understood heuristically as follows [l7|- Decay 
is due to the spreading of the wave function away from 
the initial metastable state. However, the wave packet 
keeps spreading even after tunneling out of the initial 
domain. If the evolution proceeds without any wave- 
function collapse due to measurements, the tail of the 
outgoing wave packet will reach back to the initial state. 
Since the spreading is a power-law process (the root mean 
square uncertainty in position grows proportionally to 
time), there will be a power-law tail of the wave packet 
that overlaps with the initial domain. Hence, the prob- 
ability of remaining in the initial state has a power-law 
late-time asymptotic. These considerations apply to tun- 
neling processes in field theory as well because tunneling 
occurs essentially along a one-dimensional path in field 
space, corresponding to the instanton solution. 

On the technical level, a necessary condition for the 
existence of the subexponential asymptotic is that the 
Hamiltonian of the system must be bounded (either from 
below or from above) . An elementary consideration is as 
follows. The probability of remaining in the metastable 
state is 



say by E — Eq. Using the spectral decomposition, 

/>OC 

H = / EP E dE, 

Je 



(77) 



P(t)= «>|e im |V> 



(76) 



where H is the total Hamiltonian of the system. Let us 
assume that the spectrum of H is bounded from below, 



where Pe is an orthogonal projector onto the subspace 
of energy E, we find 

{il>\j ttt W) = l e iEt (iP\P E \iP)dE= e iEt p(E)dE, 

J Eq J — OO 

where, by definition, the function p{E) identically van- 
ishes for E < Eq. Because of the nonanalyticity of p(E) 
at E = Eq, the Fourier transform of p(E) necessarily 
has a power-law asymptotic oc t~ d at t — » oo, where the 
power d is determined by the order of the (upper) nonzero 
derivative of p{E) at E = Eq. 

I conclude with some general comments regarding the 
plausibility of the age-dependent decay in cosmological 
landscape scenarios. The subexponential asymptotic was 
obtained by a quantum-mechanical consideration with- 
out regard for gravitational effects. However, gravitation 
plays a central role in vacuum decay [To| . Since the as- 
sumption of a bounded Hamiltonian is important, while 
the Hamiltonian for General Relativity is unbounded, it 
is not immediately clear that the subexponential late- 
time decay will be manifest also when the effects of grav- 
ity is taken into account. 

Another relevant consideration is the influence of 
measurements and decoherence on the vacuum decay. 
The power-law asymptotic of the survival probability 
holds only if the evolution of the wave function of the 
metastable system is unitary and proceeds according to 
the Schrdinger equation. The power-law decay can occur 
only if no wave function collapse takes place during that 
evolution. Therefore, a direct observation of the power- 
law decay is possible only if the metastable system as well 
as any decay products are perfectly isolated and do not 
have any possibility of interacting with any environment 
at least until times t ~ T. It is clear that such a perfect 
and long-lasting isolation is impossible in practice. Any 
realistic metastable system and its decay products will 
interact with an environment long before the crossover 
time T. After an interaction, the wave function will ef- 
fectively collapse back to the initial metastable state, and 
the effects of the slow spreading of the wave packet will 
be removed. 

However, one needs to be careful when applying 
quantum-mechanical considerations in the cosmological 
context. Since the potential observers of vacuum decay 
are inside the decaying field configuration, it is unclear 
whether they are able to effect a collapse of the wave 
function of the entire Hubble patch around them. Sev- 
eral points of view are possible. One could assume that 
a "measurement" of the field in the false vacuum state 
already occurs if sufficiently many gravitationally inter- 
acting macroscopic bodies are present. In that case, the 
wave function of the decaying field is continuously col- 
lapsing back to the false vacuum configuration, and so 
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it would appear that all vacuum decay is entirely inhib- 
ited due to the quantum Zeno effect (QZE), whereby a 
metastable system does not collapse when continuously 
measured. This conclusion appears implausible. On the 
other hand, it is hard to implement a measurement of 
the field values on cosmological super-horizon scales by 
any causal system. Hence, one could assume that "mea- 
surements" are absent until a tunneling event is com- 
pleted and a causally autonomous Hubble-size bubble of 
true vacuum is formed. Then one finds that the late- 
time decay asymptotic is indeed relevant to describing 
the landscape dynamics. Alternatively, one can suppose 
that "measurements" due to gravitationally induced deco- 
herence are effectively "performed" only on super-Hubble 
time and distance scales, as is the case in the decoherence 
of primordial quantum fluctuations in an inflationary uni- 
verse H [H, |2l SH. In this case, the QZE sets in only 



if the Hubble time is smaller than the time scale of onset 
of the exponential decay law. In principle, the QZE time 
scale can be estimated in a particular model of vacuum 
decay. 

Presently, I merely summarized possible viewpoints on 
the relevance of decoherence, the quantum Zeno effect, 
and subexponential decay to cosmological evolution of 
false vacuum. More work is needed to clarify this funda- 
mental issue. 
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